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USE  OF  CHANNELS  WITH  POROUS  WALLS  FOR  STUDYING  FLOWS  WHICH  OCCUR  DURING 

COMBUSTION  OF  SOLID  PROPELLANTS 


Y.  I.  Yagodkir.,  Academy  of  Sciences  USSR,  Moscow 

1 

■A. 

The  laws  of  gas  flows  in  a  channel  of  burning  charge  of  solid  pro¬ 
pellant  is  of  great  interest  for  internal  ballistics.  The  rate  of  burn¬ 
ing  of  fuel  depends  not  only  on  its  physical-chemical  properties,  but 
also  cn  the  conditions  of  gas  flew  -  laminar  or  turbulent,  and  also  on 
the  characteristics  of  turbulence  near  the  fuel’s  surface.  An  increase 
in  the  speed  of  the  gas  along  the  burning  surface  above  seme  threshold 
value  leads  to  an  increase  in  the  speed  of  burning  of  fuel.  This  effect, 
called  erosion,  is  connected,  in  the  majority  of  cases,  with  an  increase 
in  the  convective  flow  of  heat  to  the  fuels  surface.  We  propose  that  in 
charges  of  simple  form,  the  presence  of  a  threshold  speed  is  explained 
by  the  phenomenon  of  transition  from  a  laminar  mode  of  flow  to  turbulent. 
Although,  in  view  of  experimental  difficulties,  this  transition  has  not 
been  studied  under  burning-charge  conditions,  this  proposition  was  re¬ 
futed  on  the  basis  of  the  fact  that  the  value  of  threshold  Reynolds  num¬ 
ber  is  10-100  times  greater  than  the  Reynolds  number  for  transition  in 
tubes.  In  the  solution  to  this  problem,  much  aid  can  be  rendered  by  si¬ 
mulating  the  flow  and  combustion  which  occur  in  the  charge  channel  by 
using  the  channels  'with  porous  walls. 

Experimental  works  devoted  to  the  study  of  flows  in  channels  with 
porous  walls  are  very  rare.  In  the  works  of  J.  Taylor  [1]  and  V.  Veydzh- 
man  and  F.  Gevar  [2]  they  acquired  support  for  the  generality  of  dis¬ 
tribution  of  speed  in  tubes  with  porous  walls  with  one  closed  end  in 


connection  with  the  theoretical  solution  of  A.  Berman  [3]  for  a  laminar 
flow  of  viscous  liquid.  The  theoretical  distribution  of  the  axial  compo¬ 
nent  of  speed  has  the  fcrm:  < 


where  a  -  radius  of  the  tube,  v  -  speed  on  its  axis  (v~x  with  blowing 
of  gas  through  the  wail  with  a  constant  speed  vQ).  This  relationship 
is  valid  for  high  values  of  Reynolds  blowing  number  ReQ=avn/v  and  is 
a  good  approximation  with  Re0>100. 

Turbulent  flows  were  studied  by  R.  Olson  and  Ye.  Ekkert  [4]  with  the 
feed  of  gas  not  only  through  the  wail  but  also  through  one  end  of  the 
tube.  It  was  established  that  this  flow  strongly  differs  from  a  flow 
formed  only  by  blowing  through  the  wail. 
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In  the  tresent  work  we  have  continued  the  experiments  ’which  were 

^  ^  '  "*  pv»  v>  3_v*-“  4  3  1  ^  s  ]  p  0  •"PV-*  0  v»  0  ^  r»r*  pv  g  0  3.  S  *”0  J  v~ 

transition  from  a  laminar  mode  of  flow  to  turbulent  and  the  characteris¬ 
tics  of  turbulence  in  porous  tubes  with  the  feed  of  gas  through  the  wall 
Some  tests  were  dene  also  with  suction  and  blowing  through  the  leading 
end  of  the  tube  which  is  closed  by  a  porous  head.  We  used  tubes  made 
from  various  materials  (graphite,  metal,  ceramics),  with  diameters  from 
30  to  90  mm  and  with  a  length  of  up  to  1  m,  with  different  thicknesses 
of  the  walls  and  with  dimensions  of  the  pores  from  10  to  60  urn.  Measure¬ 
ments  were  done  with  the  aid  of  a  thermoanemometer  with  a  constant  tem¬ 
perature  of  the  filament,  which  had  a  diameter  of  12  pm  and  a  length  of 
1.5  mm.  The  sound  of  the  thermoanemometer  was  introduced  into  the  flow 
through  the  tube  with  an  internal  diameter  of  2.6  mm,  glued  to  the  wall 
of  the  porous  tube  (Figure  1).  Tests  were  done  with  a  pressure  within 
the  tube  from  0.2  to  1  atm. 

In  these  experiments  we  acquired  the  following  results: 

1)  With  an  increase  in  consumption  of  gas  through  the  wall  of  the 
tube,  i.e.  -with  an  increase  in  the  value  of  the  Reynolds  blowing  number 
Re„,  greater  than  some  value,  with  a  fixed  distance  x  the  sound  from  the 
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closed  end  of  the  tube,  there  occurred  in  the  flow  fluctuations  of  speed 
directly  in  its  entire  cross  section.  The  maximum  amplitude  of  the  fluc¬ 
tuations  was  found  at  a  distance  of  approximately  0.2  radius  of  the  tube 
from  the  wall.  In  tubes  with  sufficiently  fine  and  uniformly  distributed 
pores,  these  fluctuations  first  had  a  periodic  character  (Figure  2)1. 
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This  circumstance  positively  showed  the  fact  that  the  transition  was  cor. 
r.ected  with  the  hydrodynamic  stability  of  the  flow.  Figure  3  constructs 
the  dependence  of  the  mean-square  value  of  pulsations  of  the  longitudina 
component  of  speed  on  the  Reynold's  blot/ing  number.  The  graph  notes  the 
zones  of  laminar  flow  (I),  periodic  fluctuations  (II),  intermittance  of 
the  periodic  and  turbulent  fluctuations  (III)  and  turbulent  pulsations 
(IV). 

The  distribution  along  the  radius  of  the  tube  of  the  mean-square  and 
instantaneous  values  of  amplitude  of  periodic  fluctuations  of  speed  are 
shown  in  Figures  4  and  4a.  These  fluctuations  occur  in  some  relatively 


narrow  layer  which,  in  connection  with  the  theory  of  hydrodynamic  sta¬ 
bility,  can  be  called  critical.  The  frequency  of  the  pre-transition 


3 


periodic  fluctuations  changed  with  a  change  in  the  length  of  the  tube 
and  pressure  so  that  the  Strukhal  number  Sh0=af/Vg  remained  constant 
and  equalled  11  with  a  change  in  the  Reynold's  blowing  number  from  100 
so  300. 
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Figure  5.  Dependence  of  critical  Reynold's  number  on  tube's  length. 


2).  We  determined  the  critical  value  of  the  Reynold's  number  of  the 
axial  flow  Re=av/v ,  which  corresponds  to  the  beginning  of  the  occurrence 
of  speed  fluctuations.  Measurements  which  were  conducted  in  different 
tubes  showed  that  this  value  depends  on  the  dimensionless  coordinate  of 
transition  point  x/a.  The  critical  value  of  the  Reynold's  number  of  flow 
grew  together  with  a  growth  in  the  coordinate  and  reached  a  value  of 
-20,300  with  x/a=5t  (Figure  5 ) • 

If  we  accept,  according  to  A.  5erman's  theory,  a  general  distribu¬ 
tion  of  speeds  along  the  tube's  radius,  then  from  the  equation  of  con¬ 
sumption  we  can  acquire  a  dependence  of  the  Reynold's  flow  number  on  the 


coordinate  x  and  the  Reynold's  blowing  number: 

Re  =  n/?e„  — 
a 

with  the  condition  that  ReQ>100.  Using  this  dependence,  we  can  acquire 
a  coordinate  of  the  beginning  of  transition  xcr_/a  as  a  .function  of  only 
the  Reynold's  blowing  number  (Figure  6).  Let  us  note  that  this  conclu¬ 
sion  also  follows  from  the  theory  of  dimensionality. 
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Figure  5.  Dependence  of  point  of  beginning  of  transition  on  the  Reynold 

blowing  number. 

A  strong  influence  of  the  Reynold's  blowing  number  on  the  stability 
of  flow  is  determined,  apparently,  by  the  presence  of  the  transverse 
component  of  flow  speed,  which  is  not  considered  in  the  theory  of  sta¬ 
bility  of  parallel  flows.  At  small  lengths  we  also  see  as  possible  the 
influence  of  vortical  zones  near  the  channel's  head.  Recently,  V.  N. 
Varapayev,  upon  our  request,  accomplished  an  estimate  of  the  flow  in  th 
first  part  of  a  planar  channel  with  porous  walls.  With  ReQ>50  alongside 
the  head  of  the  channel,  there  are  actually  formed  two  stationary  vor¬ 
tices  (Figure  7).  They  influence  flow  on  a  length  equal,  approximately, 
to  the  width  of  the  channel  (Figure  3).  At  a  greater  length,  a  general 
distribution  of  speeds  is  established. 
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sure  5.  Distribution  of  steeds  in  first  tart  of  the  olanar  channel. 


3)  After  transition  to  a  turbulent  mode,  the  distribution  along  the 
dius  of  the  mean  sceed  remained  near  the  distribution  in  the  laminar 
ow  (Figure  9).  It  is  known  that  in  tubes  with  impermeable  walls,  ther 
curs  during  transition  a  sudden  filling  of  the  Poiseuille  distributic 
speed.  This  difference  attests  to  a  different  character  of  distribu- 
cn  of  ruisatior.  steed  and  Reynold's  stresses  in  channels  with  pcrcus 
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o -a  porous  tube. 


ana  impermeable  walls. 


The  mean-square  value  of  pulsations  of  the  longitudinal  component 
of  speed  in  the  turbulent  mode  proved  to  be  substantially  lower  than  in 
tubes  with  impermeable  ‘walls  with  the  same  values  of  Reynold's  number 
and  relative  length  of  the  tube.  For  a  comparison.  Figure  9  also  presents 
the  results  of  measurements  of  J.  Laufer  [7].  We  should  point  out  the 
particularly  large  difference  of  values  of  pulsation  speeds  near  the 
axis  and  wall  of  the  tube. 

41  With  the  feed  of  gas  through  the  front  end  of  the  tube  which  is 
closed  by  a  porous  head,  it  was  shown  that  the  transition  was  accomplishe 
with  lower  values  of  the  Reynold's  blowing  number  Reg  (Fig.  10).  Converse 
ly ,  with  suction  through  the  head  the  Reynold's  blowing  number,  which 
corresponds  to  the  beginning  of  the  transition,  increased  1.5  times  in 
the  tube  with  a  relative  length  of  x/a=11.4.  The  basic  reasonf  for  this 
growth  can  be  considered  the  removal  of  vortical  zones  in  the  tube's  head 
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In  the  second  part  of  the  work,  we  studied  the  distribution  of  the 
front  of  a  flame  in  a  uniform  propane-air  mixture  fed  within  the  tube 
through  the  porous  wall  (Fig.  1).  In  the  tests  we  determined  the  position 
and  form  of  the  flame  front  with  the  aid  of  ionization  sounding  and  pho¬ 
tographing  through  a  quartz  window  which  served  as  the  channel  head. 
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3  blowing  number  Through  the  wall 
;r  of  flow  Through  the  wall's  head. 


We  found  two  wave  modes  of  propogation  of  the  flame:  with  tangential 
waves  in  the  form  of  longitudinal  grooves  (Figure  11a)  and  with  longitu¬ 
dinal  waves  in  the  form  of  rings  (Figure  lib).  The  tangential  'waves  whic 
were  observed  at  great  distances  from  the  flame  front  to  the  wall  were 
staoior.ary,  and  for  mixtures  with  a  richer  composition  -  running  along 
The  perimeter  of  the  flame.  The  longitudinal  waves  proved  to  be  station¬ 
ary  ,  since  an  oscillogram  of  the  ionization  current  I  showed  periodic 
fluctuations,  and  the  angular  lines  on  the  flame  were  immobile.  Pictures 
of  this  form  of  fluctuations  in  Figure  lib  were  obtained  in  a  tube  made 
of  graphite  with  a  diameter  of  90  mm  with  a  low  content  in  fuel  mixture. 
The  frequency  of  fluctuations  was  f=60  Hz,  and  wavelenth  -  A  =  2  cm. 

In  this  place,  along  the  length  of  the  tube  where  the  longitudinal 
fluctuations  are  found,  the  heat  transfer  to  the  wall  increased.  The 
transition  to  a  turbulent  mode,  determined  by  the  oscillogram  :f  the 
ionization  flow,  occured  with  a  Reynold's  flow  number  of  Revi=  «mv 
coefficient  of  kinematic  viscosity  of  combustion  products)  is  less  "he 
further  the  front  of  The  flame  is  positioned  from  the  tube's  wall.  In 
This  case,  the  flow  is  less  stable  than  isothermic. 

limilar  phenomena  cf  instability  of  flame,  probably,  are  possible 
with  combustion  of  solid  fuels  in  narrow  channels.  Here,  the  speed  cf 
erosion  combustion  can  become  anemally  larger  if  the  width  cf  the  reacti 
zone  in  the  gas  phase  is  more  than  C.2  of  the  radius  of  the  channel. 


With  an  approximation  of  flame  to  wall,  it  takes  a  cylindrical  form 
and  remained  laminar  with  the  same  or  an 
flow  number  than  for  an  isothermic  flow. 

Instability  of  a  cylindrical  form  of 
follows  from  the  Landau- Mar k o h t e y n  theor; 
the  cooling  influence  of  the  wail,  we  can  say 
th-  flame  becomes  stable  for  all  wavelengths. 

form  of  stationary  longitudinal  waves  with  a  periodic  mode  over 
not  explained  by  this  theory.  In  this  case,  we  propose  that  the 
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this  case, 


e  flow  was  not  general  along  the  length  of  the  charge 
From  the  theoretical  work  of  Tszui  [11],  we  car.  conclude  a  general  dis¬ 
tribution  of  steed  and  temperature  in  the  charge,  the  flow  in  which  is 
formed  only  as  a  result  of  its  combustion. 

For  creating  a  theory  of  erosion  combustion,  we  must  further  study 


The  use  of  normal  laws  for  inpermeable  tubes  or  a  boundary  layer  cannot 
be  considered  justified. 
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